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PLASTIC  DEFORMATION  IN  BEAMS 
UNDER  SYMMETRIC  DYNAMIC  LOADS^ 

by 

J.  A.  Seller  and  P.  S#  Svmonds 
(Brown  University; 


Plastic  deformations  in  beams  caused  by  dynamic  loads 
distributed  over  finite  lengths  are  considered  in  this  paper* 

The  analysis  is  based  on  the  assumption  of  "plastic-rigid” 
behavior  so  that  the  results  should  be  valid  when  plastic  deform- 
ations are  sufficiently  large.  The  final  deformations  for 
distributed  loads  are  found  not  to  differ  much  qualitatively 
from  those  to  be  expected  for  concentrated  loads.  However, 
the  magnitude  of  the  deformation  decreases  rapidly  as  the 
length  of  the  loaded  area  Increases  from  zero,  for  given  load 
pulses*  Thus  the  simplifying  assumption  of  a mathematically 
concentrated  load  may  cause  a considerable  over-estimate  of  the 
final  deformation. 


1.  The  results  presented  in  this  paper  were  obtained  in  the 

course  of  research  sponsored  by  the  Office  of  Naval  Research 
under  Contract  N7onr-358lO  with  Brown  University. 
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PLASTIC  DEFORMATION  IN  BEAMS 
UNDER  SYMMETRIC  DYNAMIC  LOADS 
by 


1 


J.  A.  Seller  and  P.  S«  Symonds 
(Brown  University) 

1.  Introduction* 

In  a recent  paper  [1]  a method  was  described  for  cal- 
culating the  permanent  deformation  of  a bar  of  ductile  metal 
subjected  to  a transverse  impact  load  concentrated  at  the  mid- 
point of  the  bar.  The  basis  of  the  method,  originally  suggested 
by  Prager,  is  that  elastic  strains  are  neglected  entirely,  the 
bar  being  assumed  to  be  rigid  except  at  cross-sections  where  the 
limit,  cr  fully  plastic,  moment  M^  is  naintalned*  At  such  seoe 
tlons  it  is  assumed  that  plastic  strains  of  indefinitely  largo 
magnitude  can  occur*  Since  discontinuities  in  slope  angle  can 
take  place  at  such  sections,  while  the  bending  moment  remains 
essentially  constant  there,  these  sections  are  conveniently 
termed  "plastic  hinges".  This  concept  has  been  found  to  have 
valuable  applications  in  predicting  static  loads  which  produce 
failure  of  the  "plastic  collapse"  type  in  civil  engineering 
structures.  References  [2,  3]  may  be  consulted  for  descriptions 
of  these  applications,  and  for  details  concerning  the  signifi- 
cance and  calculation  of  the  limit  moment  magnitude. 

Whether  or  not  the  concept  of  localised  deformation  at 
plastic  hinge  sections  can  be  usefully  applied  in  problems  of 


1.  Numbers  in  square  brackets  refer  to  the  Bibliography  at  the 
end  of  the  paper. 
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dynamic  loading  depends  largely  on  how  much  the  energy  actually 
absorbed  during  plastic  deformation  under  a given  dynamic  load 
exceeds  the  maximum  elastic  strain  energy  which  could  be  stored 
in  the  structure*  A criterion  based  on  this  hypothesis  was 
presented  in  [l]*  Experimental  investigations  designed  to  test 
the  predictions  of  this  so-called  "plastic-rigid”  type  of  dynamic 
analysis  will  be  highly  desirable,  and  it  Is  hoped  that  the  re- 
sults of  the  present  paper  will  be  helpful  in  planning  and  inter- 
preting such  experimental  work* 

It  will  be  particularly  desirable,  in  connection  with 
experiments  and  with  practical  applications  in  general^  to  know 
how  ceT'tain  load  characteristics  may  be  expected  to  influence 
the  results.  In  a previous  paper  [4],  the  Influence  of  charac- 
teristics of  force  pulses  was  investigated,  the  force  being 
supposed  mathematically  concentrated  at  the  mid- section  of  the 
beam;  as  in  [1],  the  case  of  a uniform  beam  with  free  ends  was 
considered.  In  [4]  a simple  empirical  formula  was  suggested  by 
means  of  which  the  final  central  angle  of  deformation  produced 
by  a concentrated  force  could  be  estimated  from  knowledge  of  the 
total  impulse  and  peak  force  value;  the  shape  of  the  force-time 
curve  was  found  to  be  relatively  unimportant* 

The  present  paper  is  concerned  with  the  differences 
from  the  results  of  [4]  to  be  expected  when  the  load  is  dis- 
tributed over  a finite  length  of  the  beam  rather  than  concen- 
trated at  the  mid- section*  In  order  to  obtain  both  qualitative 
and  quantitative  comparisons  as  simply  as  possible,  the  only 
form  of  load  distribution  curve  used  is  that  indicated  in  Pig*  1* 


3 


The  total  load  P is  uniformly  distributed  over  a length  2e  where 
e £ Thus  the  parameter  specifying  the  distribution  of  the 
load  is  the  ratio  5 

A preliminary  discussion  of  this  problem  has  been  given  by 
Salvador!  and  DiMagglo  [5]«  They  used  a smoother  but  more  com- 
plicated load  distribution  fiinction,  the  load  intensity  being 
written  as: 


p(x) 


P 

2-t[l  - e“°(l  + |)1 


C|(l  + c ^)e‘° 


(1) 


where  P is  the  total  load,  x is  distance  from  midsection,  and  c 
is  a concentration  parameter  which  may  vary  from  zero  (corres- 
ponding to  a uniform  distribution  over  the  entire  beam),  to  in- 
finity (corresponding  to  concentration  of  the  load  at  the  mid- 
section). The  authors  carried  their  study  only  up  to  the  point 
of  showing  that  a sufficiently  "spread  out"  load,  characterized 
by  a value  of  the  parameter  c less  than  3.45,  would  apparently 
introduce  a qualitative  difference  in  the  behavior  o'!*  the  beam. 
This  apparent  difference  will  be  discussed  later  in  the  present 
paper.  With  a simpler  form  of  load  distribution,  as  used  in 
this  paper,  it  is  not  difficult  to  carry  the  investigation 
further,  and  to  obtain  quantitative  comparisons  between  plastic 
deformations  occurring  with  various  degrees  of  localization. 
These  quantitative  comparisons  and  more  complete  qualitative 
comparisons  between  final  plastic  deformations  are  the  main  pur- 
pose of  the  present  paper.  The  main  practical  result  is  that  a 
small  degree  of  spreading  of  the  load  over  a finite  length  of 
the  beam  leads  to  a considerable  reduction  of  the  total  deforma- 
tion at  the  mid section.  Qualitatively,  however,  the  present 


analysis  shows  that  the  final  shape  of  the  deformed  beam  will  not 
differ  radically  from  that  predicted  by  analyses  assuming  the 
force  to  be  concentrated  at  the  mid section* 

2*  Analysis  of  Earlier  Stages  of  Deformation. 

Taking  the  load  to  be  uniformly  distributed  over  the 
length  2e  as  shown  in  Pig.  1,  we  consider  the  successive  phases 
of  motion  and  deformation  which  can  occur.  We  suppose  the  total 
load  P Increases  steadily  up  to  a maximum  value  and  then  de- 
creases to  zero,  so  that  force  pulses  of  the  general  shape  shown 
in  Fig.  2 are  considered* 

Under  these  conditions,  as  the  force  increases  from 
zero  its  first  effect  is  to  cause  a translation  of  the  bar  paral- 
lel to  itself.  Then  all  points  of  the  bar  have  the  acceleration; 
a = , where  m is  the  mass  per  unit  length.  This  initial 

translation  is  termed  the  "first  phase"  of  the  motion*  Tho  con- 
stant acceleration  has  the  effect  of  a uniformly  distributed 
loading,  and  the  bending  moment  at  any  section  x [positive  bend- 
ing moments  defined  as  in  Fig.  (Vb)  ] is  given  by: 

X > e;  M(x)  = - a 

2 

X < e:  M(x)  = - a + -L  ie,  - 

2 2e  2 ’ 

where  x is  measured  from  the  middle  of  the  beam*  The  bonding 
moment  in  this  phase  has  a minimum  at  x = 0 given  by 

M(0)  = - + ^ + ^ , (2a) 

This  first  phase  ends  when  the  central  moment  reaches  the  value 
corresponding  to  a load  Pj  given  by 
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= -k- 
1-1 


Where,  for  convenience,  we  define  ^ ^ and 

In  the  "second  phase"  of  the  motion,  the  two  halves 
remain ‘rigid  but  rotate  with  respect  to  each  other,  as  in  Fig. 
(3a).  Let  denote  the  linear  acceleration  of  the  center,  and 


do  = 


denote  the  angular  acceleration  of  the  right  hand 


half  (positive  clockwise).  Assuming  to  remain  small  the 


equations  of  motion  are: 


I = m^(ao  - i ^Cq);  |(|  - |)  - Mq  = 


Solving  these  for  ag  and  a^,  we  obtain: 


So 


-[  p.(2  - I k)  - 6 ] 


= -%[3ji(l  - k)  - 12]  . 
mV 


This  second  phase  continues  until  the  moment  at  some 

additional  cross-section  reaches  the  magnitude  Mq, 

To  see  how  this  may  occur,  it  is  helpful  to  refer  to 

typical  diagrams  of  effective  load,  shear  force,  and  bending 

moment  in  the  right  hand  half  of  the  beam.  The  effective  load 

intensity  q,  with  sign  defined  as  in  Pig.  4(b),  can  be  regarded 

as  made  up  of  three  parts,  as  shown  in  Fig.  4(a);  these  are  the 

intensity  i = — 2 ^ of  the  given  load  distributed  over  the 

length  e,  the  uniform  load  intensity  -ma(— )=  -— 2 E corresoonding 

0 2 ^2  2 ^ ^ 

to  the  acceleration  a(i)  of  the  center  of  gravity  of  the  half 
beam,  and  finally  the  linearly  varying  effective  load  intonslty 
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corresponding  to  the  angular  acceleration  of  the  half-beam, 
whose  maximum  value  is  ma^  2 ” 2 ^"5  - k)  - 6],  q,  V,  M, 

curves  are  shown  in  Pigs.  (5a),  (5b),  and  (5c). 

Two  limiting  conditions  are  evident  from  the  diagrams 
of  Fig.  (5).  In  the  first  place,  unless  q(^)  is  positive,  the 
¥ curve  cannot  cut  the  axis  in  the  unloaded  position  x > e,  and 
hence  there  can  be  no  positive  moment  in  the  beam.  Hence,  in 
order  for  it  to  be  possible  to  develop  a plastic  hinge  at  a point 
where  the  maximum  moment  reaches  the  value  +Mq,  we  must  have  p, 
and  k related  so  that 


1 (1  - k)  - 6 

-1  >0 

(6) 

|i(l  - 1 k) 

> 6. 

(6a) 

Secondly,  at  the  left  hand  end  of  the  half  beam,  q(0)  must  re- 
main positive;  if  it  did  not  V would  have  negative  values  in  the 
neighborhood  of  x = 0,  and  hence  M(x)  would  decrease  below  -Mq 
in  this  region.  This  requires  the  following  relation  between 
pi  and  k: 

^ - I ti  (1  - k)  + 6 - t > 0 (7) 

or 

;i(3k  - 1)(1  - k)  < 12k.  (7a) 


It  will  then  be  evident  that  the  second  phase  may  be 
terminated  in  some  cases  because  there  is  a positive  moment  in 
the  interior  of  the  half-beam  (as  in  Fig.  5(c))  which  reaches 
the  value  at  a sufficiently  large  value  of  the  load;  and  in 
other  cases  because  the  moment  in  the  neighborhood  of  the  mid- 
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section  reaches  the  negative  limit  moment  value  -Mq  over  a finite 
length  of  the  beam,  as  in  Fig.  3(c). 

The  first  case,  occurrence  of  a ’’positive"  plastic 

hinge  in  the  interior,  is  possible  only  if  the  inequality  (6)  is 

satisfied.  From  6(a),  therefore,  it  is  evident  that  in  order 

for  this  new  hinge  to  form  in  the  interior  of  each  half-beam  it 

is  necessary  to  have:  _ 

k < I . (8) 

3 

The  second  case,  of  spreading  of  the  central  hinge,  can  occur 
only  if  inequality  (7)  is  not  satisfied,  and  hence  only  if 

k >i  (9) 

since  inequality  (7a)  will  be  satisfied,  regardless  of  p.,  if 
k <1/3. 

In  the  case  of  a concentrated  force,  [1,  4]  , the 
second  phase  ends  when  the  positive  bonding  moment  reaches  the 
value  Mq,  and  this  always  occurs  when  the  load  parameter  p reaches 
the  value  Pjj  = 22.89.  Tho  considerations  of  the  preceding  para- 
graph shov;  that  when  the  load  is  sufficiently  spread  out  (k  > 2/3), 
tho  formation  of  interior  positive  hinges  will  not  take  place  in 
this  manner.  This  is  the  main  qualitative  difference  caused  by 
loading  over  a finite  length  of  the  beam.  Instead  of  new  out- 
board hingos  making  their  appearance,  for  sufficiently  large 
k values  the  now  effect  is  tho  spreading  of  tho  central  hinge, 
which  in  the  earlier  part  of  the  motion  is  always  concentrated 
(according  to  the  basic  hypotheses  of  the  method)  at  a single 
section. 
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We  now  put  in  more  precise  terms  the  conditions  for 
(1)  occurrence  of  a positive  plastic  hinge  in  the  interior  of 
the  beam,  and  (2)  spreading  of  the  central  negative  plastic  hinge. 
(1)  Occurrence  of  lateral  plastic  hinge: 

For  e < X < -t  , the  bending  moment  is  given  by 

M(x)  = I (x  - |)  - Mq  - . (10) 

This  moment  will  have  a stationary  maximum  at  a point  \7hich  can 
be  found  by  solving  the  equation 

V = 4M  = £ - mxa.  + = 0.  (11) 

dx  2 o 2 

Using  Eqs.  (5)  this  yields  # 


(12) 


This  value  of  x is  put  in  Eq.  (10)  to  give  the  expression  for  the 
maximum  moment.  Then  the  value  of  p,  for  which  the  second  phase 
ends  is  found  by  setting  the  resulting  expression  equal  to  Mq. 

The  following  cubic  equation  in  and  k is  obtained; 


+ 36(121c^-l8k*7)iijj  + 423(5)t-4)|ij.j  + 3456  = 0.  (I3) 


For  any  particular  value  of  k,  the  real  root  of  Eq.  (I3)  then 
gives  the  value  of  p for  which  the  lateral  hinges  first  appear, 
on  the  assumption  that  the  kinematic  picture  of  Fig.  (3b)  is 
correct. 

(2)  The  central  hinge  must  begin  to  spread  if 

= q(0)  = 0. 
x=0 


(14) 
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The  critical  relation  is  therefore  obtained  by  replacing  the 
inequality  in  (7a)  by  an  equality,  which  yields 


‘‘ll  ° (3H 


12k 

- 1)(1  - k) 


(15) 


The  curves  of  p.  as  a function  of  k according  to  Eqs* 

(13)  and  (15)  are  plotted  in  Fig.  6.  It  is  found  that  the  posi- 

tive branch  of  the  curve  from  Eq.  (15)  crosses  that  from  Eq.  (I3) 
at  k = 0.361,  = 81.4.  We  conclude,  therefore  as  follows: 

(1)  If  k < 0.361  the  second  phase  ends  with  the  formation 
of  outboard  hinges.  The  initial  location  of  these  hinges  will 
be  found  from  Eq.  (12),  and  the  value  of  the  load  parameter  p,  at 
which  they  appear  is  the  real  root  of  Eq.  (13)*  Plastic  regions 
appear  in  the  manner  indicated  in  Fig.  (3b). 

(2)  If  k >0.361  the  termination  of  the  second  phase  will 
be  due  to  the  spreading  of  the  central  plastic  hinge  into  a 
finite  region  in  which  the  moment  has  the  constant  value  -Mq. 

The  load  parameter  at  which  this  effect  begins  is  given  by  Eq. 

(14) ,  and  the  corresponding  type  of  deformation  is  indicated  in 
Fig.  (3c). 

This  critical  value  of  k = 0. 36I  corresponds  to  the 
critical  value  c = 3*^5  found  by  Salvador!  and  DiMaggio  for  the 
smooth  load  distribution  curve  assumed  by  them,  Eq.  (1).  However, 
they  refer  to  "splitting  of  the  central  hinge"  as  the  effect  re- 
sulting when  Eq.  (l4)  is  satisfied  before  the  lateral  hinges 
form.  They  state  that  when  Eq.  (l4)  is  satisfied,  "the  moment  M 
becomes  minimum ,at  x = 0 and  the  central  hinge  splits  into  two 
symmetrical  hinges  which  wander  outward,  plasticizing  a central 
portion  of  the  beam".  This  suggests  that  two  distinct  hinges  are 
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present  in  place  of  the  original  central  hinge,  with  the  portion 
of  the  beam  between  them  subjected  to  a moment  less  in  magnitude 
than  Mor  and  hence  moving  as  a rigid  body.  This  could  only  occur, 
however,  if  the  load  intensity  has  a relative  minimum  at  the  cen- 
ter, which  is  not  the  case  in  the  distribution  curve  assumed  by 
them.  In  the  usual  type  of  load  distribution  met  with  in  impact 
problems  in  which  the  load  intensity  is  either  constant  on  a 
relative  maximum  at  the  midpoint,  the  effect  is  not  a "splitting" 
but  a spreading  of  the  central  hinge.  In  other  words,  unless 
the  load  Intensity  increases  (near  x = 0)  with  distance  from  the 
midpoint,  the  effect  is  the  development  of  a "hinge"  of  finite 
length,  with  the  moment  constant  at  the  value  -Mq  in  a finite 
region  which  increases  in  size  as*  the  load  increases. 

For  the  analysis  of  the  motion  after  the  close  of  the 
second  phase  we  assume  first  that  the  load  distribution  parameter 
k ^ e/^  is  less  than  0. 36I  and  the  1 that  k exceeds  0.361. 

3«  Analysis  for  k < 0.l6l. 

In  this  case  the  lateral  hinges  appear  while  the  cen- 
tral hinge  is  still  localized  at  the  center  section.  The  appro- 
priate free  body  diagrams  for  the  right  hand  half  are. shown,  in 
Fig.  (7),  the  distance  of  the  lateral  hinge  from  the  center  being 
denoted  by  Xj^.  It  is  assumed  that  Xj^  > e.  The  equations  of 
motion,  as  in  [l,  4]  are 

I = "«h(ao  - j 

I ( ^ - f ) - 2M„  = ^ (l6b) 

0 = mCi  - Xh)[aj^  Xh^“l^ 


(17a) 
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(17b) 


0 

Mq  — ni(^  ” ®1* 

The  fifth  equation  relating  the  five  unknowns  Eq,  ai,  a^, 
and  is  that  which  expresses  the  discontinuity  of  linear  accel- 
eration at  the  hinge  location  Xj^,  which  is  in  general  a function 
of  time.  This  is  (See  [1,  4]) 

dx. 

®o  - *h®o  “ t®l  “l^'^  * *h^^=  dt^^“o  “ 


where  «q,  are  the  angular  velocities  of  the  segments  to  the 
left  and  to  the  right,  respectively,  of  the  hinge,  i. e. , 


a 


0 


dwQ 

“ dt  ’ 


(19) 


Equations  (l6)  - (l8)  yield  the  following  system  of  equations: 


dH)  . 
®o  “ dt 


3M, 


- k)  - 8] 


m 


12Mr 


m- 


1^3(1  - 


(20a) 

(20b) 


= — !1S-5  [^(45  - 3k)  - 24] 
2mlh^ 

6Mq 

~ m^2(i  - 


(20c) 

( 2oa) 


where 


K * 


(20e) 

(21) 


These  equations  hold  so  long  as  the  bending  moment  has 
a magnitude  less  than  Mq  everywhere  except  at  the  sections  x = 0 
and  X = x^.  The  effective  load,  shear,  and  moment  diagrams  for 
the  segments  to  the  left  and  right  of  the  moving  hinge.  Fig.  (8) 
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indicate  that  the  critical  condition  will  be  a spreading  of  the 
central  hinge.  This  will  occur  when  the  following  condition  is 


satisfied 


dx 


x=0 


= q(0)  = 0. 


(22) 


x=0 


From  Eq.  (16)  this  leads  to  the  condition 


^ “ (^-k)(3k  - O * 

The  use  of  the  limiting  condition  requires  that  Eqs.  ( 2d)  first 
be  solved,  yielding  (among  other  quantities),  as  a function 
of  p,;  these  solutions  will  be  valid  as  long  as  n is  less  than 
the  value  computed  from  Eq.  ( 2]i). 

If  any  given  finite  force  pulse  is  replaced  by  a rec- 
tangular pulse  of  the  same  peak  force  and  the  same  area,  i.e. , 
total  impulse,  the  final  deformations  computed  for  the  rectangular 
pulse  will  exceed  the  true  ones  for  the  original  pulse,  but  for 
many  purposes  the  results  will  be  useful  estimates  (e.g. , in 
error  by  15  or  20  percent)^  see  [4],  Therefore  to  obtain  quan- 
titative comparisons  of  total  deformations  with  various  values 
of  k we  have  assumed  a rectangular  force  pulse  shape,  Pig.  (2b). 
Since  the  details  of  the  method  of  analysis  exactly  parallel 
those  in  [4],  only  the  final  results  are  shown  in  Fig.  (12)  where 
curves  of  the  deformation  parameters  are  plotted  as 

functions  of  = PjjjVMq  for  a series  of  values  of  the  distribu- 
tion parameter  k.  These  results  and  others  are  discussed  further 
in  tho  final  section  of  the  paper. 
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4-t  Analysis  for  k > Q« 

In  this  case  the  central  hinge  begins  to  spread  before 
lateral  hinges  appear.  At  time  t the  central  hinge  will  have 
spread  out  a distance  d from  the  center;  in  the  region  -d  < x < d 
the  moment  is  constant  at  the  value  The  equations  of  motion 

for  the  segment  of  length  i - d are  obtained  from  the  free  body 
diagrams  of  Fig.  (9)»  and  are  found  to  be; 


■^(e  - d)  = m( 't  • d)  [a^  ('f  -d)ci^] 


(24a) 


i(e  - d)C^  -.2^]  - d)3a’  (24b) 


2e 


where  a^  is  the  acceleration  at  x = d and  is  the  angular  accel- 
eration of  the  segment  of  length  -t  - d. 

Equations  (24)  yield  expressions  for  the  accelerations 


as  follows; 


^0  = 


<*0  = 


Mo 


m^(l  - K) 

12Mq 

m^^d  - X)^ 


5 ■ W<4  . X - 3k)  - 6] 

[^(k  - k)(l  - k)  - I] 


Where 


(25a) 

(25b) 

(250 


To  determine  d,  we  use  the  fact  that  the  central  hinge  spreads 
at  a rate  Just  sufficient  to  prevent  the  moment  near  x = d from 
decreasing  below  the  value  -M©.  The  condition  is  that 


dx^ 


x=d 


= q(d) 


(26) 


Using  Eq.  (2va)  and  simplifying,  it  is  found  that  X is  the  follow- 
ing function  of  ^ and  k; 
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X - nQk  - D(1  - k)  - 12k  ^ (27) 

2^(1  - k) 

Plots  of  vs.  X for  various  values  of  k are  given  in 
Pig.  (10).  Thus  for  k > 0. 36I  the  critical  portion  of  the  beam 
(-  X ^ X-t  £ x)  has  the  constant  moment  -Mq,  so  that  the  accelera- 
tions of  all  elements  are  proportional  to  the  local  load  intensity 
P/2e.  The  two  segments  of  length  -t-ds-td-X)  rotate  as  if 
rigid  bars  hinged  at  the  sections  x = + d,  Fig,  (3c).  So  long 
as  the  moment  in  the  outer  segments  rer^ins  less  than  Mq  in  mag- 
nitude, this  is  the  correct  picture.  The  load  magnitude  at 
which  this  picture  ceases  to  be  the  correct  one  is  that  at  which 
new  "positive"  hinges  form  in  the  Interior  of  the  outer  segments. 
The  maximum  bending  moment  in  these  segments  can  be  investigated 
in  the  same  manner  as  was  done  for  the  case  of  a localized  cen- 
tral hinge.  The  following  equation  was  obtained  by  setting  the 
maximum  bending  moment  in  one  segment  d ^ x ^ ^ equal  to  Mq.  It 
is  analogous  to  Eq.  (16)  but  depends  now  upon  X = d/^  in  addition 
to  k. 


2 - - X)^(l  - X)^[ti(k  - X)(2  X - 3k)-12k]  ^ ^(k  - X)^ 

4k[3^(k  - X)(l  - k)  - 12k]^ 

+ ti^(k  - X)^(l  - X)3  ^ 

6k[3ji(k  - X)(l  - k)  - 12k]3 


To  determine  the  load  at  \i?hlch  the  additional  hinges  first  appear, 
the  above  equation  must  be  solved  simultaneously  with  Eq,  (27), 
Alternatively,  for  any  chosen  value  of  k ( >0.361)  plots  of  p.  vs. 

X from  Eq,  (28)  and  Eq,  (27)  can  be  drawn,  and  the  intersection 
of  the  two  curves  yields  the  initial  value  of  p,  in  question.  As 
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one  example,  the  case  k = 0*5  has  been  worked  out*  In  this  case, 
when  the  two  equations  are  combined  we  obtain 

- 192|i^  + Il52ii2  - 110592  = 0.  (29) 

This  has  solution  (x  = 186,  corresponding  to  X = 0*242;  the  new 
hinge  appears  at  x = 0*660<t* 

5*  Numerical  Results  and  Conclusions. 

The  type  of  load  distribution  considered  here,  with 
uniform  load  intensity  over  part  of  the  beam  and  zero  load  over 
the  rest,  is  Intended  less  as  a realistic  loading  than  as  the 
simplest  type  of  distributed  load  to  handle  analytically*  It 
seems  intuitively  clear  that  the  qualitative  results  deduced  from 
this  simple  case  will  apply  to  general  cases  of  loads  spread  over 
finite  areas.  The  main  qualitative  result  is,  in  fact,  that  the 
general  picture  of  the  final  deformation  is  exactly  the  same  as 
that  in  the  case  of  concentrated  loading  unless  the  length  over 
which  the  load  is  applied  exceeds  about  a third  of  the  total 
length  of  the  beam*  Only  when  k exceeds  the  value  0. 361  does  a 
qualitative  difference  appear*  This  consists  of  a spreading  of 
the  central  hinge,  so  that  instead  of  a plastic  hinge  localized 
at  the  midsection  there  is  a finite  length  of  the  beam  in  which 
the  bending  moment  has  the  magnitude  Mq*  However  in  real  beams 
a hinge  cannot  be  localized  at  a single  cross-section  in  any 
case  since  this  would  imply  Infinite  strains;  the  plastic  region 
is  always  spread  over  a finite  region  even  when  the  load  is 
strongly  localized*  Thus  the  final  plastic  deformation  of  a beam 
subjected  to  a sufficiently  large  distributed  load  will  not  in 
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any  case  differ  radically  from  that  predicted  by  an  analysis 
using  a mathematically  concentrated  load. 

As  far  as  quantitative  comparisons  are  concerned  it  is 
of  particular  Interest  to  Investigate  how  a small  degree  of  load 
spreading  modifies  the  magnitude  of  the  final  deformation  pre- 
dicted by  the  simpler  analysis  for  a concentrated  force,  since 
in  practice  loads  are  always  applied  over  finite  areas. 

To  provide  information  of  this  kind  calculations  have 

been  made  of  the  final  central  angle  Oq  produced  by  a series  of 

* 

loads  with  k values  ranging  from  sero  to  0.30.  The  force-time 
curve  was  taken,  for  maximum  simplicity,  to  have  a square-wave 
shape,  as  in  Fig,  2(b).  For  this  type  of  force  pulse  the  lateral 
hinges  move  Instantaneously  to  positions  which  they  maintain  so 
long  as  the  load  remains  constant.  The  lateral  hinge  coordinate 
is  found  as  the  solution  of  the  following  equation: 


^ ^ (1  - 0^ 


(30) 


This  is  Eq.  (20e)  with  the  right  hand  side  set  equal  to  zero. 
Figure  11  shows  curves  of  ^ plotted  against  ^ for  values  of  k 
ranging  from  zero  up  to  O.3O.  The  calculation  of  the  further 
increments  of  deformation  which  occur  after  the  load  drops  to 
zero  are  made  in  the  manner  described  in  [4]. 

Figure  12  shows  curves  of  the  dimensionless  parameters 
ni't3©o/MQT^  plotted  as  a function  of  the  dimensionless  load  para- 
meter = PjnVMoj  for  k equal  to  0,  0.10,  0.20,  and  O.3O.  Figure 
13  shows  curves  of  m'?^©o/MQT^  plotted  against  k for  various  values 
of  M'jj.  These  curves  show  that  when  tho  maximum  load  is  large 
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(e.g*y  > *+0)  the  deformation  falls  off  very  rapidly.  For 
example,  with  = 40,  the  deformation  parameter  drops  from  about 
800  for  a concentrated  load  (k  = 0)  to  about  600  for  k = 0*0^* 
This  apparently  large  reduction  caused  by  the  distribution  of  the 
load  over  a finite  length  may  be  somewhat  exaggerated,  as  com- 
pared with  practical  cases,  because  of  the  special  conditions  of 
the  present  calculation,  i*e«,  the  assumption  of  a rectangular 
force  pulse  and  of  a constant  load  intensity  over  the  loaded  seg- 
ment of  the  beam.  However  the  results  indicate  that  the  assump- 
tion of  a mathematically  concentrated  load  may  lead  to  an  appre- 
ciable overestimate  of  the  deformation,  since  the  actual  load  is 
spread  over  a small  but  finite  length  of  the  beam#  This  result 
should  be  of  significance  in  evaluating  experimental  results. 

It  should  perhaps  also  be  noted  that  strain-hardening  will  cause 
a further  reduction  in  the  permanent  deformation;  this  effect 
will  presumably  be  more  important  the  higher  the  degree  of  con- 
centration of  the  load. 
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Pig.  4(b) 

Iffactlve  load  components 
in  the  second  phase 
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Pig.  5(c) 


Typical  diagrams  of  effective 
load  Intensity  and  associated 
shear  and  bending  moment  dia- 
grams In  second  phase. 
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Fr»e  body  diagram  of  sagmont  to 
right  of  eantral  plastic  region 
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Diagrams  of  eff active  load,  shear 
and  bending  moment  in  the  third 
I^se,  with  concentrated  central 
hinge 
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Fig,  10  Dcpendmo*  of  %rldth  d of  contral  hlngo  oa  loodl  P for. 

rarloas  valuos  of  lo«d  distribution  parsnster  k - s/v 


Plg»  12  Final  eantral  angle  of 
deformation  caused  by 
various  values  of  the 
load  distribution 
parameter  k.  as  func- 
tion of  maximum  load  P 
of  a rectangular 
force  pulse.  j 
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